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PREFACE 


This  report  documents  the  work  performed 
by  the  Communications  Croup  of  The  University  of 
Alabama  in  Huntsville  for  the  U.S.  Army  Missile 
Command  on  pulse  compression  signals  for  tracking 
applications,  under  Contract  DAAHOl-7 6-C-0991 . 

Work  continues  on  the  other  tasks  specified 
In  the  contract  and  this  research  will  be  documented 
in  subsequent  reports. 
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ABSTRACT 


This  report  reviews  pulse  compression  signals  and  their 
applicability  to  tracking  radars.  The  emphasis  is  on  binary  phase 
coded  signals,  such  as  Barker  codes  and  linear  maximal  length 
sequences.  The  properties  of  complementary  sequences  are  reviewed 
and  techniques  for  utilizing  these  sequences  in  tracking  radars 
are  discussed.  Linear  frequency  modulated  signals,  polyphase  codes 
and  signals  with  simultaneous  amplitude  and  phase  modulation  are 
also  considered. 
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1.  Introduction 


Pulse  compression  Involves  the  transmission  of  a long  coded  pulse 
and  the  processing  of  the  received  echo  to  obtain  a relatively  narrow 
pulse.  This  technique  combines  the  advantages  of  high  range  resolution 
as  In  a narrow  pulse  and  of  good  doppler  resolution  and  high  average 
power  as  In  a wide  pulse.  Pulse  compression  Is  usually  achieved  by  the 
use  of  a matched  filter  whose  output.  In  the  absence  of  doppler  shift 
and  noise,  is  the  autocorrelation  function  of  the  transmitted  signal. 
Ideally,  the  autocorrelation  roust  be  Iropulse-like  with  a large  value 
at  zero  lag  and  zero  elsewhere,  and  the  degradation  due  to  doppler  shift 
must  be  small.  The  presence  of  secondary  peaks,  or  side  lobes,  severely 
limit  the  performance  of  the  radar.  In  the  case  of  tracking  radars,  how- 
ever, the  compressed  pulse  may  be  allowed  large  side  lobes  provided  they 
fall  outside  the  range  bin  being  tracked.  A pulse  compression  signal  Is 
a pulse  of  finite  width  with  a high  bandwidth- time  product.  For  a given 
pulse  width,  the  bandwidth  may  be  Increased  by  phase  and/or  amplitude 
modulation. 

The  goal  of  the  task  reported  here  Is  to  review  and  analyze  the 
literature  In  view  of  selecting  a pulse  compression  signal  for  tracking 
applications  of  a pulsed  radar.  The  performance  criteria  are:  (1)  the 

peak-to-side  lobe  ratio  of  the  autocorrelation  function,  (11)  the  dis- 
tance between  the  main  peak  and  the  first  large  secondary  peak,  and 
(111)  the  degradation  In  the  compressed  pulse  due  to  small  doppler 


mismatches . 


The  review  of  phase  modulated  signals  Includes  continuous  phase 
modulation  and  discrete  phase  modulation.  Linear  and  nonlinear  frequency 
modulations  constitute  continuous  phase  modulation.  Discrete  phase  modu- 
lation consists  of  polyphase  codes  and  binary  phase  codes.  Barker  codes, 
shift  register  generated  sequences,  and  complementary  series  are  the 
most  Important  binary  phase  codes  that  are  considered.  Huffman  codes 
and  generalized  Barker  codes  Involve  both  amplitude  and  phase  modulation. 
Side  lobe  reduction  may  also  be  affected  by  processing  the  matched  fil- 
ter output  through  a weighting  filter.  Such  filters,  suitable  for  cer- 
tain types  of  signals,  are  also  discussed.  The  final  section  of  the 
report  contains  a comparative  analysis  of  the  results  available  In  the 


literature 
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‘ 2.  Signal  Modulation  and  Pulse  Compression  [5.  10.  14.  471 

i 

In  a pulsed  radar,  the  transmitted  signal  Is  of  the  general  form, 
s(t;a),)  » rect^(t)  a(t)  cos  (u^t  - <Kt))  (2.1) 

where 

a(t)  Is  the  amplitude  modulation, 

(|)(t)  Is  the  phase  modulation, 
o)j  - 2irf,, 

fg  Is  the  carrier  frequency, 
rect^(t)  » 1 for  |tl  < y , 

* 0,  otherwise, 

and  T Is  the  width  of  the  pulse. 

The  return  signal  from  a target  at  range  R and  moving  with  a radial 
velocity  V (range  rate)  Is 

r(t)  - k rect^(t  - t)  a(t  - t)  cos  { (Wj  - ‘*J^)(t  - t)  - 4)(t  - t)} 

- k s(t  - t;  Ug  - w^)  (2.2) 

2R 

where  T • ~ la  the  delay, 

2vWj 

• — - — Is  the  doppler  shift, 

k Is  an  attenuation  factor, 
c Is  the  velocity  of  propagation, 
and  v Is  positive  for  a target  moving  away  from  the  radar. 

Equation  (2.2)  Is  derived  using  a first  order  approximation  to  the 
round-trip  delay  due  to  a varying  target  range. 


The  received  signal  is  contaminated  by  noise  generated  in  the 
atmosphere  or  in  the  antenna  and  the  receiver.  Therefore,  some  form 
of  signal  processing  is  essential  prior  to  detection  or  target  parame- 
ter estimation.  Assuming  white  noise,  the  matched  filter  is  the 
optimal  linear  filter  which  maximizes  the  signal-to-nolse  ratio  at  its 
output  [14,  18,  57].  The  impulse  response  of  such  a filter  is 

h(t)  - s*(-t;  U),  - 0)^  (2.3) 

where  is  the  selected  doppler  match  and  is  generally  equal  to  zero. 

When  the  doppler  shift  of  the  received  signal  is  o)^  , the  output  of  the 
matched  filter  y(t)  is  given  by 

y(t)  - k4>  (t  - T)  + (|)  (t  - T)  (2.4) 

89  ilS 

where  <|)  (•)  is  the  autocorrelation  of  s(t),  and  <j)  (•)  is  the  cross- 
es ns 

correlation  between  n(t)  and  s(t).  Assuming  that  the  signal  and  the 
noise  are  sufficiently  uncorrelated,  y(t)  will  exhibit  a peak  at  t=T, 
which  is  the  desired  delay  and  is  a measure  of  the  target  range.  The 
peak  value  is  proportional  to  the  energy  in  the  signal  and  the  width 
of  the  processed  pulse  is  inversely  proportional  to  the  bandwidth  of 
the  signal.  A doppler  mismatch  alters  the  output  shape  and  a two- 
dimensional  correlation  function  must  be  Introduced.  Such  a function 
is  the  ambiguity  function  which  will  be  discussed  In  the  next  section. 

A good  range  resolution  requires  that  the  processed  signal  be  a 
narrow  spike,  which  can  be  achieved  by  using  a signal  of  large  bandwidth. 

One  of  the  means  of  obtaining  a large  bandwidth  is  the  reduction  in  the 
duration  of  the  pulse.  However,  since  the  signal  detectability  and 
accuracy  of  measurements  depend  on  the  energy  of  the  signal,  the  transmitter 
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power  must  be  raised  to  keep  the  energy  constant.  Peak  power  limitations 
and  efficient  usage  of  the  average  power  set  a lower  limit  on  the  pulse 
. duration.  The  problem,  therefore,  Is  to  increase  the  bandwidth  of  the 

i signal  Independent  of  the  pulse  duration.  A monochromatic  RF  pulse  Is 

defined  as  one  whose  amplitude  and  phase  modulations,  a(t)  and  4>(t),  are 
both  constant.  It  can  be  easily  shown  that  the  bandwidth  of  such  a 
signal  Is 

B - Y ’ (2.5) 

Thus  the  bandwidth  can  be  Increased  only  by  reducing  the  pulse  duration,  ! 

I 

and  the  tlme-bandwldth  (BT)  product  Is  equal  to  one.  In  a pulse  com- 
pression radar,  the  signal  Is  designed  to  have  a high  BT  product  where 
a large  bandwidth  is  achieved  without  sacrificing  the  pulse  duration 
[4,  7,  31,  47,  52].  This  is  accomplished  by  phase  modulation,  amplitude 
modulation,  or  both.  In  such  a situation,  the  matched  filter  converts  a 
wide  pulse  at  the  Input  to  a narrow  pulse  at  the  output  and  Is  called 
the  pulse  compression  filter.  Since  the  width  of  the  input  pulse  is  T 

and  that  of  the  output  Is  , the  ratio  of  the  input  pulse  width  to  the 

B 

output  pulse  width  Is  BT.  This  ratio  Is  also  known  as  the  compression 
ratio. 

The  most  common  method  for  Increasing  the  BT  product  Is  phase 
modulation  (PM).  The  variation  of  <|)(t)  may  be  continuous  as  In  frequency 
modulation  (linear  or  nonlinear  FM)  or  discrete  (phase-coded  signals). 

] 


3.  Ambiguity  Function 


The  two  most  important  functions  of  a radar  are  detection  of  a 
target  and  the  estimation  of  its  parameters,  such  as  position  and  ve- 
locity. Target  detection  is  usually  performed  by  subjecting  the  matched 
filter  (MF)  output  to  a thresholding  device.  When  this  output  exceeds 
the  threshold,  a target  is  detected  and  the  time  at  which  crossing 
occurs  is  a measure  of  the  target  range.  Because  of  the  presence  of 
noise,  detection  may  occur  even  in  the  absence  of  a target.  Therefore, 
the  threshold  must  be  selected  to  satisfy  the  specified  probability  of 
false  alarm.  Then,  the  probability  of  detection  can  be  increased  by 
increasing  the  peak  value  of  the  MF  output  which,  as  mentioned  earlier, 
is  proportional  to  the  energy  in  the  pulse. 

For  a given  peak,  the  accuracy  of  estimation  of  the  parameters 
depends  upon  the  shape  of  the  MF  output  [47 ] . Assuming  no  noise  inter- 
ference and  zero  doppler  mismatch,  the  output  is  the  autocorrelation 
function  of  the  transmitted  pulse.  Figure  3.1  shows  the  autocorrelation 
functions  of  three  types  of  radar  signals,  where  the  last  two  include 
side  lobes  in  addition  to  the  main  lobe.  If  the  MF  output  in  Figure 
3.1a  is  passed  through  a thresholding  device,  the  signal  exceeds  the 
threshold  for  a large  period  of  time.  Since  the  target  range  is  esti- 
mated from  the  time  instant  at  which  the  threshold  is  crossed,  the 
spread  of  the  response  introduces  an  uncertainty  or  ambiguity  as  to  the 
exact  target  range.  The  degree  of  ambiguity  is  smaller  if  the  main 
lobe  is  narrow.  Further,  if  there  are  two  targets  whose  range  separa- 
tion corresponds  to  a delay  which  is  less  than  the  width  of  the  main 
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lobe,  it  Is  Impossible  to  resolve  the  two  targets.  Figure  3.1b  shows 
a response  where  the  main  lobe  is  narrow  enough  for  the  desired  sepa- 
rability. However,  because  of  the  high  side  lobes  which  may  also  ex- 
ceed the  threshold,  the  target  may  be  at  any  one  of  several  ranges. 

It  may  also  not  be  possible  to  distinguish  the  return  due  to  a single 
target  from  a superposition  of  returns  due  to  several  targets.  An  ideal 
waveform  would  be  one  whose  autocorrelation  function  has  one  high  narrow 
peak  and  is  zero  elsewhere.  Such  waveforms  cannot  be  realized  in 
practice  and,  besides,  the  presence  of  noise  rules  out  the  feasibility 
of  an  ideal  response.  A more  practical  response  is  depicted  in  Figure 
3.1c  which  features  a narrow  peak  and  a long  tail  of  side  lobes.  Here 
again,  a weak  target  in  the  vicinity  of  a strong  target  will  be  masked 
out  by  the  side  lobes.  Therefore,  for  a given  value  for  the  peak,  the 
desired  autocorrelation  should  have  a narrow  main  lobe  and  small  side 
lobes . 

Thus  far,  the  discussion  assumes  that  there  was  no  doppler  mis- 
match. However,  since  the  target  range  rate  is  not  known  exactly,  the 
matching  signal  will  have  a different  doppler  shift  compared  to  the 
target  return.  This  could  result  not  only  in  the  reduction  of  Che  peak 
value  but  also  in  a time  shift,  thus  promoting  inaccuracies  in  the 
measurements.  Therefore,  a two-dimensional  autocorrelation  function 
called  ambiguity  function  must  be  defined  for  waveform  evaluation  in 
terms  of  ambiguity,  resolution,  and  doppler  mismatch. 

Before  defining  the  ambiguity  function,  it  is  useful  to  consider 
Che  complex  representation  of  radar  signals.  The  transmitted  signal 
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s(t)  given  In  Equation  (2.1)  can  be  written  aa  [5,  10>  14,  51] 

8(t;o),)  - Re{iKt)}  (3.1) 

where  Re  {*}  denotes  "real  part  of"  and  i])(t)  is  a complex  signal  given 
by 

jWjt 

i|;(t)  - u(t)  e (3.2) 

and  u(t)  - a(t)  (3.3) 

Similarly,  assuming  zero  delay,  the  target  return  r(t)  can  be  written  as 


r(t)  - Re{A(t)}  (3.4) 

jWot 

A(t)  - v(t)  e (3.5) 

.»i  v(0  . k.(t) 

-jw.t 

- ku(t)  e . (3.7) 


The  analysis  could  be  carried  out  either  with  the  real  signals  s(t)  and 
r(t)  or  equivalently,  with  the  corresponding  complex  signals  u(t)  and 
v(t).  Assuming  that  the  filter  Is  matched  to  the  signal  with  zero 
doppler,  the  Impulse  response  of  the  filter  Is  [57] 

h(t)  - u*(-t)  . (3.8) 

The  output  of  the  matched  filter  is 

w(t)  ■ v(t)  * 1i(t) 

-jw.t 

■ k u(t)  u*(t  - t)  e dt  . (3.9) 

A two-dimensional  correlation  function  X(t,u)^)  Is  Introduced,  defined 


as 
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jo)dt 


X(T,u)d)  - J u(t)  u*  (t  - T)  e “ dt 
The  filter  response  is  easily  seen  to  be 
w(T)  - kX(T,-u>d) 


(3.10) 


(3.11) 


and  the  real  response  of  the  matched  filter  is 

Re{kX(T,-o)d)  e ® } . 

The  two-dimensional  function  |X(x ,0)^)1  is  known  as  the  uncertainty 
function  and  |x(T,u)d)|*  is  called  the  ambiguity  function  [47,  64].  The 
ambiguity  function  shows  the  performance  of  the  signal  with  respect  to 
ambiguity,  resolution  and  the  effect  of  doppler  shift.  Throughout  this 
analysis,  it  was  assumed  that  the  matching  signal  had  zero  doppler. 
However,  if  the  filter  is  matched  to  a doppler  of  , then  the  variable 
o>d  in  the  ambiguity  function  represents  the  difference  between  0)^  and 
the  actual  doppler  shift.  Some  of  the  Important  properties  of  the 
ambiguity  function  are  given  below  without  proof  [5,  47]. 

(a)  If  U(u))  is  the  Fourier  transform  of  u(t),  then  an  alternate 

expression  for  X(t,o).)  would  be 

a 

X(T,0)d)  / U*((D)U(u>-Wd)  e^'^'^du)  . (3.12) 

(b)  The  ambiguity  surface  is  symmetrical  with  respect  to  the 
origin;  l.e., 

|X(-T,  -0Jd)r*  [X(T, 01^)1^  . (3.13) 

(c)  The  value  at  the  origin  Is  proportional  to  the  square  of  the 
signal  energy. 
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|X(0,0)1*  - |/|u(t)|"dt]' 

■ 4 (signal  energy)^  (3.14) 

(d)  The  absolute  maximum  value  of  the  ambiguity  surface  is  at  the 
origin;  thus 

|X(T,Wj)|*  i 1X(0,0)P  . (3.15) 

(e)  The  total  volume  under  the  ambiguity  surface  is 

//|X(T,u)^)|^dx  du)^  = 2h1x(0,0)12  . (3.16) 

The  above  equation  states  that  the  total  volume  is  a function  only  of 

the  signal  energy,  and  whatever  Is  done  to  make  |X(t,o)^)|^  small  in  one 

part  of  the  T,u.  plane  will  result  in  an  increase  someifhere  else.  Thus, 
a 

if  one  chooses  a signal  whose  ambiguity  has  a narrow  central  spike,  the 
bulk  of  the  volume  must  appear  as  side  lobes,  possibly  masking  targets 
«dilch  are  relatively  far  away.  If  It  Is  Important  to  avoid  such  masking, 
the  volume  must  be  concentrated  about  the  origin,  resulting  In  a broad 
central  peak  with  Its  Inherent  disadvantages. 

The  matched  filter  receiver  Is  optimized  for  signal  detection  In 
white  noise,  but  may  have  the  disadvantage  of  relatively  high  side  lobes. 
In  this  case,  a weighting  filter  Is  used  In  conjunction  with  the  matched 
filter  to  reduce  the  side  lobes  (14,  51].  Such  a deliberate  mismatch 
Inevitably  degrades  noise  performance,  but  It  Is  conceivable  that  a 
good  compromise  can  be  found  between  improvement  In  resolution  and  de- 


terioration In  detection  performance.  The  most  common  reason  for  mis- 
matching Is  the  suppression  of  range  side  lobes.  Specifically,  the 


heights  of  the  side  lobes  are  reduced  at  the  cost  of  widening  and 
reducing  the  main  lobe. 

In  the  application  being  considered  in  this  study,  the  radar  is 
assumed  to  be  in  the  tracking  mode  and,  therefore,  good  estimates  of 
the  range  and  the  range  rate  are  available.  Thus  only  the  matched  fil- 
ter output  corresponding  to  a few  range  cells  around  the  estimated 
range  are  considered,  the  rest  of  it  being  essentially  Ignored.  Large 
side  lobes  can  therefore  be  allowed  to  exist  outside  this  Interval  of 
interest.  This  means  that  only  a small  portion  of  the  ambiguity  sur- 
face around  the  origin  needs  to  be  examined  while  evaluating  the 
performance  of  a given  waveform. 


A.  Linear  and  Nonlinear  FM 


The  Instantaneous  frequency  of  the  complex  envelope  u(t)  Is  defined 


as 


(4.1) 


PM  waveforms  for  which  the  phase  4>(t)  Is  nonlinear  but  where  It  Is  con- 
venient to  express  the  modulation  of  the  Instantaneous  frequency  are 
called  FM  waveforms.  Linear  FM  [32,  37,  44]  Is  defined  by  an  Instan- 
taneous frequency  which  Is  a linear  function  of  time. 


fl(t)  - ^ wt  for  f . t . I 
and  a corresponding  phase  modulation. 


<Kt) 


The  complex  modulation  function  u(t)  Is  given  by 

u(t)  . ^ 


(4.2) 


(4.3) 


(4.4) 


B B 

The  frequency  sweep  of  this  signal  Is  - ~2  to  — , where  the  bandwidth 


B Is 


B 


UT 

2ti 


(4.5) 


It  can  be  easily  shown  that  the  uncertainty  function  of  this  waveform 
Is  [5,  14,  20] 


|X(T,0).)  I 


sln| 

(MT  - U)  .> 

U 

|t 

fi-M' 

L TJ 

( 2 

lit 

-“d 

2 

(4.6) 


and  Is  depicted  In  Figure  4.1.  For  small  values  of  T,  Ignoring  the 
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relatively  slowly  varying  factor  ^1  - the  function  1X(t,u^)|  can 


be  approximated  as 


|X(T,0)^)|  S T 1 Sa  (WT  - W^)>1 
where  Sa(x) 


d' ' ~ * I 2 
alnx 


(4.7) 


Thus,  any  cut  parallel  to  the  delay  or  doppler  axis  yields  a 
|Sa(x)|  function,  with  the  peak  occurring  at 


MT  - u)^  • 0 


or 


(4.8) 


In  other  tiords,  a doppler  mismatch  of  results  In  an  equivalent  time 
“d 

delay  of  , thus  reducing  the  accuracy  of  the  range  measurement.  The 

factor  ^1  - decreases  the  amplitude  as  one  moves  away  from  the 

doppler  axis.  Thus  a doppler  mismatch  results  In  a reduction  of  the 
peak  and  delays  It  by  an  amount  proportional  to  the  mismatch.  In  the 
vicinity  of  the  origin,  the  response  along  the  delay  axis  can  be  approxi- 
mated as 


. Sin— CPll/2) 
2 / (yTT/2) 


(yTx/2) 

sin  ttBt 
hBt 


(4.9) 


The  half  power  width  of  the  central  peak  Is  of  the  order  of  (1/B) . Since 
the  Input  pulse  width  is  T,  the  compression  factor  is  BT,  the  tlme- 
bandwldth  product.  Similarly,  It  can  be  shown  that  the  half  power  width 
In  doppler  la  about  (1/T) . The  autocorrelation  structure  described  by 
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(4,9)  also  Indicates  Che  existence  of  severe  side  lobes.  In  fact,  the 
height  of  Che  first  side  lobe  Is  within  -13.5  dB  of  the  central  peak. 

Two  drawbacks  of  linear  FM  are:  (1)  the  range-doppler  coupling 

as  suggested  by  (4.8),  and  (2)  the  existence  of  relatively  large  range 
side  lobes.  In  many  applications,  a good  estimate  of  the  doppler  Is 
available,  and  the  range  estimate  can  be  corrected  to  account  for 
doppler  mismatch.  Otherwise,  It  Is  necessary  to  use  an  entire  bank  of 
filters,  each  matched  to  a particular  value  of  doppler.  This,  of 
course.  Is  also  necessary  if  a doppler  measurement  Is  contemplated. 

The  range  side  lobes  can  be  suppressed  using  weighting  filters 
[14,  55,  63],  thus  creating  an  Intentional  mismatch.  As  mentioned 
earlier,  weighting  filters  achieve  side  lobe  suppreslon  at  the  cost  of 
reduced  resolution  and  signal  detectability.  The  loss  In  resolution 
Is  a result  of  broadening  of  the  main  lobe  and  the  reduction  In  the 
central  peak  (hence,  the  slgnal-to-nolse  ratio)  causes  the  loss  In 
detectability.  The  Dolph-Chebyshev  weighting  [14]  is  optimum  In  the 
sense  that  It  has  the  smallest  main  lobe  expansion  for  a specified  side 
lobe  suppression.  However,  It  Is  not  physically  realizable  and  several 
approximations  to  It  have  been  suggested.  The  most  Important  among 
these  are  the  Taylor  approximation,  Hamming  weighting,  and  the  general 
cosine-power  weighting.  A sixth  order  Taylor  weighting  designed  for  a 
peak  side  lobe  level  of  -40  dB  also  results  In  a main  lobe  widening  by 
a factor  of  1.41  and  In  a loss  In  the  slgnal-to-nolse  ratio  of  1.2dB. 

The  Hamming  weighting  is  considerably  less  difficult  to  realize  than  the 
Taylor  weighting  and  Is  a particular  case  of  a cosine-squared  function  on 


i 
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a baseline  pedestal.  The  performance  figures  for  Hamming  weighting  are: 
peak  side  lobe  level  of  -42.8  dB,  main  lobe  widening  factor  of  1.47, 
and  loss  In  slgnal-to-nolse  ratio  of  1.34  dB.  Reduced  side  lobes  can 
also  be  achieved  by  amplitude  weighting  or  amplitude  modulation  of  the 
transmitted  signal. 

The  strong  coupling  between  range  and  range  rate  can  be  attributed 
to  substantial  match  In  the  Instantaneous  frequencies  of  a delayed, 
doppler-shlfted  signal  and  that  of  the  matching  signal  as  seen  In 
Figure  4.2. 


Figure  4.2  Instantaneous  Frequency  of  Linear  FM  Signal 
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One  possible  method  of  lessening  this  match  is  the  use  of  nonlinear 
I modulation  such  as  quadratic  FM  and  piecewise  linear  FM  [47].  Even 

) though  the  peaks  due  to  range-doppler  coupling  are  reduced,  high  ridges 


appear  elsewhere  because  the  volume  under  the  ambiguity  surface  Is 
constant. 
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5.  Polyphase  Codes  [21,  27,  40,  41.  51] 


The  previous  section  was  concerned  with  continuous  frequency  (or 
phase)  modulation.  In  this  and  the  following  sections,  a class  of  sig- 
nals called  the  phase  coded  signals  will  be  discussed.  These  waveforms 
feature  ordered  sequences  that  are  Impressed,  at  fixed  time  Intervals, 
on  the  phase  of  the  carrier.  In  other  words,  the  phase  modulation  <J>(t) 
Is  assumed  to  be  piecewise  constant. 


<l>(t)  “ (J)jj  for  (n  - 1)6  i t i n6  , n = 1,2,  , N (5.1) 


where  T » N6  . 


(5.2) 


Polyphase  codes  with  zero  periodic  correlation  except  at  zero  shift 
were  developed  by  Helmlller  [27].  These  codes  are  of  length  N = p^  where 
p Is  a prime  number  and  p different  phases  corresponding  to  the  p*"^  roots 
of  unity  are  required  to  generate  the  codes.  The  restriction  that  p be 
prime  was  later  removed  by  Frank  and  Zadoff  [22].  Frank  [21]  described 
a method  for  generating  non-perlodlc  polyphase  codes  using  the  same  con- 
struction as  for  periodic  codes.  The  code  of  length  N = p^  is  described 
by  a p X p matrix: 


(p-1) 
2(p-l) 
3 (p-1) 


(5.3) 


0 (p-1)  (2(p-l)  (p-1)' 


20 


The  elements  represent  multiplying  coefficients  of  a basic  phase  angle 
(2TTk/p)  where  k Is  an  Integer  and  Is  relatively  prime  to  p (k  Is  usually 
equal  to  1) . The  non-perlodlc  sequence  Is  obtained  by  placing  the  rows 
side  by  side,  yielding  a sequence  containing  p^  elements.  For  example, 
let  p ■ 4;  then  the  matrix  Is 

0 0 0 0 

0 12  3 

(5.4) 

0 2 4 6 

0 3 6 9 

and  the  basic  phase  angle  Is 


The  sequence  of  phases  (<|»^,  n»l,2,  ••••,  16}  can  be  written  as 

{0,  0,  0,  0,  0,  e,  20,  30,  0,  20,  40,  60,  0,  30,  60,  90}  . 

As  mentioned  before,  the  codes  were  originally  designed  to  give 
zero  periodic  correlation,  and  this  optimum  side  lobe  performance  Is  not 
obtained  with  non-perlodlc  correlation.  Cyclic  shifting  of  the  codes 
given  above  results  In  new  codes  with  different  non-perlodlc  autocorre- 
lation structures.  But  It  has  been  found  that  the  sequence  as  given 
above  possesses  the  most  desirable  structure. 

The  phase  variation  of  the  Frank  polyphase  code  can  be  considered 
as  a piecewise  constant  approximation  to  a parabola;  thus  the  Frank 
polyphase  code  Is  an  approximation  to  linear  FM. 

Frank  has  conjectured  that  the  maximum  side  lobe  Is  the  vector  sum 
of  p/2  (for  p even)  or  (p  + l)/2  (for  p odd)  unit  vectors  where  the 
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vectors  are  separated  by  2ir/p.  This  implies  that  for  large  p,  the  peak 
side  lobe  amplitude  approaches  p/ir  and  the  peak-to-slde  lobe  ratio 
approaches  irp  [14].  The  actual  ratios  for  values  of  p from  3 to  8 are 
given  in  Table  5.1. 

TABLE  5.1  Feak-to-Slde  Lobe  Ratios  for  Polyphase  Codes 


p 

N 

Peak  Side  Lobe 
Amplitude 

Peak-to-Slde  Lobe 
Ratio  (dB) 

3 

9 

1.0 

19.08 

4 

16 

1.4 

21.16 

5 

25 

1.6 

23.88 

6 

36 

2.0 

25.11 

7 

49 

2.25 

26.76 

8 

64 

2.6 

27.82 

i 


Matched  filter  outputs  for  doppler  shifted  polyphase  codes  of 
length  N ■■  100  are  shown  In  Figures  5.1(a)  - 5.1(g)  for  different 
values  of  doppler  shifts  [14].  It  can  be  seen  from  these  figures  that 
the  polyphase  codes  have  doppler  degradations  somewhat  similar  to  those 
of  linear  FM.  The  main  peak  shifts  Its  position  as  a function  of 
doppler  shift  and  its  value  decreases  as  the  doppler  shift  Increases. 

In  addition,  there  are  secondary  peaks  which  are  cyclic  as  functions 
of  doppler  shift.  At  higher  dopplers,  an  Image  signal  also  develops 
so  that  the  output  breaks  up  Into  two  similar  signals. 


Figure  5.1  Matched  Filter  Outputs  for  Polyphase  Codes 
with  Different  Doppler  Mismatches 
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Figure  5.1  Matched  Filter  Outputs  for  Polyphase  Codes 
with  Different  Doppler  Mismatches 


Other  three-phase  codes  consisting  only  of  the  cube  roots  of  unity 
were  Investigated  by  DeLong  [16].  Existence  of  such  sequences  has  not  been 
established  for  lengths  greater  than  9.  Four-phase  codes  have  also  been 
considered  because  of  their  simplicity  of  Implementation.  Four-phase  approxi- 
mations to  Frank  codes  of  length  greater  than  16  and  to  linear  FM  have  been 
studied. 

Sequences  constructed  from  two  or  more  phases  and  possessing  the 
property  that  the  magnitude  of  the  autocorrelation  at  any  non-zero  delay  Is 
less  than  one  were  Investigated  by  Golomb  and  Scholtz  [26].  Such  sequences 
are  called  Generalized  Barker  Sequences,  and  a four-phase  sequence  of  length 
IS  has  been  described.  Also  Included  In  this  paper  Is  an  exhaustive  tabu- 
lation of  six-phase  sequences  of  lengths  N ■ 3 to  N > 12.  Existence  of 
longer  sequences  has  not  been  established  even  though  there  appear  to  be 
no  theoretical  constraints  which  would  preclude  them. 

Codes  which  Involve  amplitude  modulation  In  addition  to  phase  modu- 
lation have  also  been  Investigated.  Some  of  these  codes,  such  as  Huffman 
codes  and  ternary  sequences,  will  be  described  In  Section  9.  When 
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the  number  of  distinct  phases  Involved  In  a polyphase  code  Is  2,  the 

I 

! code  reduces  to  a binary  phase  code.  Because  of  the  extensive  research 

performed  on  binary  codes,  they  are  treated  separately  In  the  following 
sections. 
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6.  Barker  Codes 


The  simplest  form  of  phase  coding  Is  the  binary  or  phase-reversal 
code  in  which  the  phase  (j>^  is  either  zero  or  180°.  The  complex  repre- 
sentation of  such  a signal  is  given  by. 


ip(t)  = u(t)  e 


(6.1) 


where 


u(t)  = e for  (n-l)6  < t S n6 


0 or  180” 


n = 1,2, 


(6.2) 


(6.3) 


and  T 


(6.4) 


Combination  of  (6.2)  and  (6.3)  yields 


y a R < 


(6.5) 


where  R (t)  = 1 for  (n-l)6  < t < n6 
n 


= 0 elsewhere 


(6.6) 


and  a = ± 1 . 

n 


(6.7) 


It  can  be  easily  shown  that  the  autocorrelation  of  such  a waveform  is 
made  up  of  triangular  pulses,  and  hence  X(T,0)  may  be  obtained  by  cal- 
culating the  function  X(k6,0)  for  integral  values  of  k and  combining 
these  points  by  straight  lines.  The  discrete  samples  of  the  auto- 


correlation are  given  by 


X(k6,0) 


■f ' V 


n+lk| 


(6.8) 
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The  Barker  codes  or  Barker  sequences  are  a family  of  binary  se- 
quences characterized  by  [5,  14,  15,  43,  51,  59] 

X(k6,0)  “ N for  k = 0 

= 0,  ± 1 for  k 0 • (6.9) 

Because  of  this  desirable  property,  they  are  also  called  perfect  words 
and  optimal  sequences.  The  size  of  this  family  of  codes  is  restricted; 
in  fact.  Barker  sequences  are  known  to  exist  only  for  seven  lengths. 

These  are  given  in  Table  6.1,  where  + and  - represent  +1  and  -1, 

respectively. 


It  has  been  proved  that  Barker  sequences  do  not  exist  for  odd  lengths 
greater  chan  13  and  that  the  existence  of  even  length  Barker  code  re- 
quires that  N be  a perfect  square  (39,  53,  58,  59).  Further,  it  has 
been  shown  that  no  even  length  Barker  sequences  exist  for  N < 6084  [58]. 
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Figure  6.1  shows  the  magnitudes  of  the  matched  filter  outputs  for 
different  values  of  doppler  mismatch.  The  figure  shows  slices  at  doppler 
Intervals  of  (2ir/4T),  with  the  zero  doppler  slice  In  the  front.  As  can 
be  easily  seen,  the  response  deteriorates  rapidly  for  large  values  of 
doppler  mismatch.  Because  of  the  high  peaks  exhibited  at  non-zero 
dopplers,  these  codes  are  not  desirable  even  if  a bank  of  matched  fil- 
ters is  used  for  all  possible  dopplers. 

Even  for  the  longest  available  Barker  code  (N  = 13),  the  peak  side 
lobe  at  zero  doppler  mismatch  Is  only  22.3  dB  below  the  main  peak.  In 
a multi-target  environment,  such  high  side  lobes  are  undesirable.  As 
In  the  case  of  linear  FM,  one  way  to  suppress  the  side  looes  Is  In  the  use 
of  a weighting  filter.  The  geometric  similarity  between  the  main  lobe 
and  the  side  lobes  (they  are  triangles  of  equal  width)  is  used  to  syn- 
thesize a weighting  filter  using  delay  lines.  Figure  6.2(a)  shows  the 
autocorrelation  function  p(t)  of  the  Barker  sequence  of  length  13.  Each 
of  the  side  lobes  can  be  independently  reduced  to  any  desired  level  by 
adding  weighted  and  delayed  replicas  of  the  function  itself  [20,  38]. 
Consider  a 2K  1 order  network  with  impulse  response, 

K 

h(t)  - ^ 3^6(t  - 2k6)  (6.10) 

k-  -K 

where  6(')  is  the  Dirac  delta  function.  If  the  input  to  this  filter  is 
P(t),  then  it  can  be  easily  shown  that  its  output  Is  given  by 
K 

- 53  ’ 

k--  K 


g(t) 


(6.11) 


t 


Figure  6.2  Barker  Code  Autocorrelation  Before  and 
After  Weighting 
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It  is  desirable  to  maintain  the  evenness  of  the  response;  thus, 

■ ^-k  • 

It  Is  now  required  to  choose  K + 1 weighting  factors  3j , 
to  provide  the  desired  side  lobe  suppression.  One  possible  criterion 
Is  to  reduce  the  first  K side  lobes  to  zero  while  maintaining  a value 
of  13  for  the  main  lobe.  For  example,  if  N =13  and  K=6,  this  cri- 
terion results  in  7 linear  equations  in  the  seven  unknown  wleghtlng 
factors.  Solution  of  these  equations  yields  the  following  values: 

= 1.047722  , 3,  = -0.040733  , = -0.045572  , 3j  = -0.050094 

3^  = -0.054269  , 3j  = -0.058066  , 3g  = -0.061461  . 

The  weighting  filter  can  be  easily  mechanized  with  a delay  line  network 
as  shown  in  Figure  6.3.  The  output  of  the  filter  g(t)  is  shown  in 
Figure  6.2(b).  As  can  be  seen,  the  first  six  side  lobes  have  been  com- 
pletely suppressed  but  the  weighting  filter  has  introduced  additional 
side  lobes  at  delays  larger  than  136.  The  main  lobe-to-side  lobe  ratio 
of  this  signal  is  about  32.4  dB,  an  improvement  of  10  dB  over  the  non- 
weighted  signal.  More  important,  the  side  lobes  at  smaller  delays 
are  completely  cancelled,  making  it  attractive  for  tracking  applications. 
A 25th  order  network  designed  to  suppress  the  first  12  side  lobes  yields 
a main  lobe-to-side  lobe  ratio  of  45  dB.  As  mentioned  earlier,  the  use 
of  a weighting  filter  destroys  the  optimality  of  the  matched  filter, 
resulting  in  a loss  in  detection  capability  (signal-to-noise  ratio).  It 
has  been  shown  that  the  maximum  loss  for  the  Barker  code  of  length  13  is 
0.25  dB. 

Weighting  filters  can  also  be  designed  in  the  frequency  domain  [46] . 

The  Fourier  transform  of  the  autocorrelation  function  P(t)  can  be  shown  to  be 


33 


P(co)  - 6 Sa*(w6/2)  [l2  + • (6.13) 

The  transform  of  the  desired  portion  (the  main  lobe)  Is 

G(o))  - 136  Sa^ (0)6/2)  . (6.14) 

Thus  the  desired  signal  can  be  extracted  by  Inverse  filtering;  l.e., 
by  passing  P(t)  through  a filter  whose  transfer  function  Is 

H(u))  - G(o))/P(o)) 

" , sin  13~o)6  ' (6.15) 

sin  0)6 

Weighting  filters  have  been  designed  to  approximate  this  transfer 
function.  Delay  line  filters  based  on  this  approximation  can  be  obtained 
from  a truncated  Fourier  series  expansion  of  H(o)) . Other  approximations 
Including  the  use  of  feedback  structures  have  also  been  Investigated. 

Although  one  can  suppress  the  side  lobes  of  Barker  codes,  these 
codes  allow  only  low  pulse  compression  factors  because  of  the  limited 
length.  Longer  codes  can  be  achieved  by  phase  coding  within  each  seg- 
ment of  one  Barker  code  with  another  Barker  code.  Such  sequences  have 
been  called  combined  Barker  sequences,  multistage  Barker  codes,  or 
Barker  squared  codes  [29,  30].  Each  pulse  of  the  new  code  may  In  turn 
be  replaced  by  an  entire  code,  and  so  forth.  By  using  relatively  short 
codes  of  lengths  Nj,  Nj,  Nj,  •••*,  one  obtains  a large  overall  code 
length  of  NjNjNj  ••••  . Figure  6.4  shows  the  autocorrelation  function 
of  the  Barker  code  of  length  13  where  each  segment  Is  coded  with  a 
Barker  code  of  length  4.  When  the  roles  of  the  codes  of  length  13  and 
4 are  reversed,  the  resulting  autocorrelation  Is  as  shown  In  Figure 
6,3.  In  either  case,  a compression  ratio  of  52  has  been  achieved,  but 
the  side  lobes  are  considerably  stronger  than  those  of  a Barker  code. 
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Figure  6.4  Autocorrelation  of  a 13-bit  Barker  Code  each  bit  of 
which  is  a 4-bit  Barker  Code 
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Figure  6.5  Autocorrelation  of  a 4-bit  Barker  Code  each  bit  of 
which  Is  a 13-bit  Barker  Code 
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Here  again,  side  lobe  suppression  can  be  effected  by  the  use  of  a 
weighting  filter.  Either  of  the  two  methods  mentioned  earlier  can  be 
used  for  this  purpose.  The  applicability  of  the  second  method  can  be 
proved  as  follows  [46].  Consider  a Barker  code  consisting  of  Nj  pulses 
of  duration  6 each.  Then,  assuming  positive  side  lobes  (i.e.,  2,5 

or  13),  the  autocorrelation  of  this  code  is 

(Ni  - l)/2 

P,(T)  = (N,  - l)d(T)  + d(T  - 2n6)  (6.16) 

n=-(Nj  -D/2 

where  d(T)  Is  the  autocorrelation  of  an  individual  pulse.  If  this 
code  Is  repeated  in  the  form  of  a positive  side  lobe  Barker  code  of 
length  Nj  , the  autocorrelation  of  the  new  code  will  be 

(Nz  - l)/2 

P2<t)  = (Nj  - 1)P,(T)  + ^ Pi(T  - 2n6)  . (6.17) 

n=  -(Nj-  l)/2 


By  taking  the  Fourier  transforms  of  Pj(t)  and  Pz(t),  it  can  be  shown 


that 


Pj(ui)  = D(a>) 


1 + 


sin  Nj  0)6 
sin  0)6 


= D(o))F,(o)) 


and  ?2(w)  = Pj(w) 


sin  N2  0)6 

N,  - 1 + — : — 1 — 

2 sin  0)6 


= P,(o))F2(o)) 

= D(0))F,(o))F2(0)) 


(6.18) 


(6.19) 

(T) 


The  desired  signal  d(x)  can  therefore  be  obtained  by  passing  P^ 
through  a filter  whose  transfer  function  is 


Approximations  to  this  transfer  function  can  be  used,  as  before,  to 
reduce  the  range  side  lobes. 


A shift  register  sequence  generator  consists  of  a shift  register 
working  In  conjunction  with  appropriate  logic,  which  feeds  back  a 
logical  combination  of  the  states  of  two  or  more  of  Its  stages  to  Its 
Input.  A shift  register  generator  (SRG)  of  degree  n consists  of  n 
storage  elements  (n  stages),  a clock  pulse  generator,  and  a feedback 
logic  circuit  as  shown  in  Figure  7.1.  When  a clock  pulse  arrives,  the 
contents  of  each  stage  are  transferred  to  the  next  stage.  The  new 
state  of  the  first  stage  Is  an  arbitrary  Boolean  function  of  the  pre- 
vious contents  of  the  register.  For  each  stage,  the  off  state  is 
symbolized  by  0 and  the  state  by  1.  The  binary  output  of  the  feed- 
back logic  unit  is  a Bo'  an  function  of  n binary  Inputs,  and  there  are 
, n 

2^  different  Boolean  functions  for  a given  n.  A shift  register  generated 
sequence  is  therefore  specified  by  the  particular  Boolean  function  chosen 


Figure  7.1  The  General  Shift  Register  Generator 


38 

and  the  initial  state  of  the  register.  It  can  be  shown  that  the  output 
sequence  of  an  n stage  generator  Is  periodic  with  period  p S 2*'  - 1 and 
that  for  every  n,  there  actually  exist  feedback  arrangements  which  lead 
to  the  maximum  period  N = 2^  - 1.  Sequences  with  the  maximum  period 
are  called  maximum  length  sequences. 

Denote  by  x^  the  state  of  the  l^*'  stage  and  let  f(x,,  Xj.  x^) 

denote  the  Boolean  feedback  function.  The  shift  register  is  called 
linear  if  this  function  can  be  expressed  as 

f(x,,  X2,  , x^)  = C,Xj  ® CjXj  © C3X3  ® ®^n*n 

where  each  constant  is  either  zero  or  one  and  where  the  symbol  @ 

denotes  addition  modulo  2.  Otherwise,  the  shift  register  is  nonlinear. 

I\  tl  2 ^ 

Since  there  are  2 ways  to  pick  the  n constants  only  2 of  the  2 
shift  registers  with  n stages  are  linear.'  An  example  of  a simple  5- 
stage  linear  SRG  is  whown  in  Figure  7.2  where 

“ Cj  * 1 and  Cj=C3=C^=0 
and 

fCXj,  Xj,  Xj,  x^,  Xj)  = X2@Xj  (7.2) 

Linear  maximum  length  (MI.)  shift  register  generated  sequences  or 
linear  m-sequences  have  been  investigated  in  great  detail.  It  has 


Figure  7.2  Five  Stage  Linear  SRG 
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been  proved  that  such  sequences  exist  for  every  register  length  n.  The 
linear  feedback  functions  which  yield  m-sequences  have  been  extensively 
tabulated  [9,  17,  65]  and  a part  of  It  Is  shown  In  Table  7.1.  The  use 
of  this  table  will  be  demonstrated  with  an  example.  Consider  a 5-stage 
generator,  n = 5 for  which  the  maximum  period  lsN=2^-l*31.  The 
table  lists  3 possible  feedback  connections: 

(1)  [5,2]  which  corresponds  to  a feedback  connection  described  by 

f(Xj,  Xj,  Xj,  x^,  Xj)  * Xj © Xj  (7.3) 

(11)  [5, 4, 3, 2]  which  corresponds  to 

f(x,,  x^,  Xj,  x^,  Xj)  = Xj @ Xj © x^ © Xj  (7.4) 

and  (111)  [5, 4, 2,1]  which  corresponds  to 

f(Xj,  X2,  Xj,  X^,  Xj)  = Xj©X2©X^©Xj  (7.5) 

Figure  7.2  shows  the  feedback  connections  for  the  first  of  these  gene- 
rators. In  addition  to  the  connections  listed  In  the  table,  there  exist 
other  combinations  which  produce  sequences  which  are  mirror  Images  of 
those  obtained  using  the  table.  For  example.  If  the  table  shows  the 
connection 

I 

[n,  t,,  tj,  tj  ] , 

the  mirror  Image  sequence  Is  obtained  from  the  connection 
[n,  n-tj,  n-tj,  n-tj,  ] . 

It  should  be  noted,  too,  that  for  a given  connection,  different  Initial 
conditions  produce  the  same  sequence  with  different  shifts.  The  output 
of  the  shift  register  generator  Is  a sequence  of  zeros  and  ones.  The 

} 

I 

1 

i 

i 

3 


Number  of 
stages 


Code 

length 

Feedback  Connections 

3 

[2.1J 

7 

[3,1] 

15 

[4.1] 

31 

[5.2]  [5. 4. 3, 2]  [5,4, 2,1] 

63 

[6,1]  [6, 5, 2,1]  [6, 5, 3, 2] 

127 

[7.1]  [7,3]  [7, 3. 2,1]  [7, 4, 3, 2] 
[7, 6, 4, 2]  [7, 6, 3.1]  [7, 6, 5, 2] 

[7. 6. 5. 4. 2.1]  [7, 5, 4, 3, 2,1] 

255 

(8, 4, 3, 2]  [8,6,5,31  [8, 6, 5, 2] 

[8. 5. 3.1]  [8, 6, 5,1]  [8, 7, 6,1] 

[8. 7. 6. 5. 2.1]  [8, 6,4, 3, 2,1] 

511 

[9.4]  [9. 6, 4, 3]  [9, 8, 5, 4]  [9,8,4. 1] 
[9.5, 3, 2]  [9, 8, 6, 5]  [9,8, 7,2] 

[9, 6. 5, 4, 2,1]  [9, 7, 6, 4, 3,1] 

[9, 8. 7. 6, 5, 3] 

1023 

[10,3]  [10,8,3.2]  [10,4,3,1]  [10,8,5,1] 
[10,8,5,4]  [10,9,4,1]  [10,8,4,3] 
[10,5,3,2]  [10,5,2,1]  [10,9,4,2] 

2047 

[11.1]  [11,8,5,2]  [11,7.3,2]  [11,5,3,2] 
[11,10,3,2]  [11,6,5,1]  [11,5,3,1] 

[11.9.4.1]  [11,8,6,2]  [11,9.8,3] 

4095 

[12.6.4.1]  [12,9,3,2]  [12,11,10,5.2,1] 

[12.11.6.4.2.1]  [12,11,9,7,6,5] 

[12.11.9.5.3.1]  [12,11,9,8,7,4] 
[12,11,9,7,6,5]  [12,9,8,3,2,1] 

[12.10.9.8.6.2] 

8191 

[13.4.3.1]  [13,10,9,7,5,4] 

[13.11.8.7.4.1]  [13,12,8,7,6,5] 

[13.9.8.7.5.1]  [13,12,6,5.4,3] 
[13,12,11,9.5,3]  [13,12,11,5,2,1] 
[13,12,9,8,4,2]  [13.8,7,4,3,2] 

16383 

[14.12.2.1]  [14,13.4,2]  [14,13,11,9] 

[14.10.6.1]  [14,11,6,1]  [14,12,11,1] 
[14,6,4,2]  [14,11,9,6,5,2] 
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binary  phase  modulated  signal  is  obtained  by  using  phase  angles  of 
zero  (a^^  ■ + 1)  and  180*  (a^  = - 1)  for  zero  and  one,  respectively. 

The  fundamental  properties  of  a binary  phase  coded  linear  m-sequence 
generated  by  an  n-stage  generator  are  [5,  12,  15,  23]: 

the  sequence  is  periodic  with  period  N = 2”  - 1, 
in  each  period,  the  difference  between  the  number  of 
ones  and  the  number  of  minus-ones  is  one, 
half  the  runs  (consecutive  states  of  same  kind)  are  of 
length  one,  one-fourth  of  length  two,  etc., 
as  a consequence  of  linearity,  if  an  m-sequence  is  added 
modulo  - 2 to  a shifted  version  of  itself,  the  resulting 
sum  sequence  is  a shifted  version  of  the  original  sequence, 

the  periodic  correlation  is  two-valued;  i.e., 

N 

S Vi+k  ° ” If  1 = 0,  t N,  ± 2N,  etc.  (7.6) 

k»  1 

“ - 1,  otherwise. 

Non-periodic  m-sequences  are  obtained  by  retaining  Just  one  period 
of  the  periodic  version.  Different  sequences  can  be  obtained  by  trun- 
cating the  output  after  N samples  and  using  different  initial  conditions. 
For  an  n-stage  register,  (2”  - 1)  different  initial  conditions  are  pos- 
sible since  the  all-zero  state  is  not  allowed.  Figures  7.3(b)  - 7.3(h) 
show  the  effectof  the  inlt  ial  condit  Ion  on  the  autocorrelation  function  of 
the  m-sequence  produced  by  the  3-stage  generator  shown  in  Figure  7.3(a) 
[14,  49].  As  can  be  easily  seen,  the  peak  side  lobe  and  its  proximity 
to  the  main  peak  vary  widely  with  the  initial  condition.  Similarly, 
sequences  generated  using  different  feedback  connections  also  have 


(i) 

(11) 

(ill) 

(Iv) 

(v) 


43 


i; 


I 

[ 

i 


t 


! 


different  autocorrelation  functions.  There  is  no  theoretical  tool 
available  to  aid  in  the  selection  of  the  feedback  connection  and  the 
initial  condition  that  will  result  in  a desirable  autocorrelation 
function.  To  study  the  effect  of  doppler  mismatch,  the  correlation 
functions  were  computed  for  an  m-sequence  generated  by  a 5-stage 
generator  with  feedback  connection  [5,3]  and  initial  condition 
(1,  0,  1,  0,  1).  Figure  7.4  shows  a cutaway  section  of  the  uncer- 
tainty function  for  doppler  limits  of 

2v  . ^ 2tt 

" 106  ‘^d  106  • 

For  purposes  of  comparison,  a few  slices  of  this  function  are  shown 

in  Figures  7.5(a)  - 7.5(g).  Only  positive  values  of  x and  oj,  are 

a 

shown  in  these  figures  since  the  uncertainty  function  is  a symmetric 
function  of  both  T and  OJ^.  For  large  values  of  N,  the  peak-to-side 
lobe  ratio  at  zero  doppler  approaches  /n  (compare  with  N for  Barker 
codes  and  ir/^  for  Frank  polyphase  codes).  As  mentioned  earlier,  the 
range  side  lobe  structure  is  different  for  different  initial  condi- 
tions. Table  7.2  shows  the  initial  conditions  for  some  of  the  linear 
ro-sequence  generators  which  yield  the  largest  peak-to-slde  lobe  ratio 
[43,  54]. 

The  range  side  lobes  may  be  further  reduced  by  the  use  of  a 
weighting  filter  similar  to  the  one  employed  In  the  case  of  Barker 
codes.  The  performance  of  m-sequences  in  the  presence  of  doppler  mis- 
match is  clearly  superior  to  that  of  the  Barker  code  since  it  has  no 
appreciable  spurious  peaks. 

For  an  n-stage  SRG,  2^  linear  feedback  combinations  are  possible. 
Only  a few  of  these  produce  maximal  length  sequences,  others  resulting 


i 


J 
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In  non-maximal  length  sequences  whose  period  p < 2*^  - 1.  For  a given 
feedback  connection,  unlike  In  the  case  of  the  m-sequences , different 
Initial  conditions  produce  sequences  of  different  periods.  For  a 
given  non-maxlmal  generator,  there  Is  a longest  period  available  and 
It  Is  obtained  by  loading  the  register  with  all  zeros  except  for 
either  the  first  or  last  stage  [5].  Because  of  this  Initial  loading, 
such  a sequence  Is  called  the  Impulse  response  sequence.  Very  little 
work  has  been  done  on  the  use  of  non-maxlmal  length  sequence  In  pulse- 
doppler  radars. 

The  generation  of  nonlinear  sequences  can  be  considered  as  falling 
Into  two  cases:  (1)  the  case  where  the  feedback  logic  Is  nonlinear, 

and  (11)  the  case  where  the  output  of  a basic  linear  generator  Is 
altered  by  nonlinear  means.  The  advantage  of  nonlinear  feedback  Is 
chat  many  more  sequences  of  maximum  period  are  possible  for  a given 
number  of  stages.  The  major  disadvantage,  however.  Is  the  lack  of  a 
suitable  theory  for  theoretical  analysis  of  nonlinear  generators.  It 
has  been  shown  that  there  are  2^^**  ' - n)  nonlinear  sequences  of 
maximum  period  for  an  n-stage  register.  In  the  second  method  of 
generating  nonlinear  sequence,  linear  m-sequences  are  combined  to 
form  composite  codes.  Even  though  these  codes  are  not  maximal  length 
sequences,  they  have  other  Interesting  properties  which  warrant  their 
Investigation.  The  Gold  codes  [17]  are  generated  by  modulo  - 2 addition 
of  a pair  of  linear  m-sequences  of  Che  same  period.  The  new  sequence 
has  the  same  period  as  Che  original  sequences  but  Is  non-maxlmal.  The 
Gold  codes  are  usful  for  generating  a large  number  of  codes  of  a given 
length.  Syncopated-reglster  generators  [17]  use  a technique  that 
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multiplexes  two  sequences  to  produce  a sequence  of  much  higher  period. 
Two  m-sequences  are  added  modulo -2,  except  that  one  of  the  generators 
has  Its  clock  shifted  by  half  a period  with  respect  to  the  other  clock. 
If  the  two  generators  have  r and  s stages  (r  ^ s) , then  the  period  of 
the  composite  code  Is  (2^  - 1)(2®  - 1).  The  JPL  ranging  codes  [17] 
are  similar  to  the  Gold  codes  but  are  constructed  by  modulo  - 2 addition 
of  two  or  more  m-sequences  whose  lengths  are  relatively  prime  to  one 
another.  The  period  of  the  composite  sequence  Is  the  product  of  the 
periods  of  the  original  sequences. 

One  other  known  class  of  binary  sequence  whose  periodic  correla- 
t^lon  Is  maximum  at  zero  delay  and  minus-one  elsewhere  Is  the  Legendre 
sequence  or  the  Perron  sequence  [24,  25].  Perron  sequences  are  of 
period  N =*  4m-l  when  this  number  is  a prime.  The  binary  phase  coded 
signal  of  this  type  Is  defined  as 


(7.7) 


where  the  Legendre  symbol  (— ) is  defined  as 

N 


1 r ^ 

■ i 


1,  if  there  is  an  Integer  x for  which  x^  = l(mod  N) 

otherwise  . (7.8) 


For  example. 


If  m = 3 and  N = 11,  it  can  be  easily  shown  that 
+1,  If  1 Is  among  0,1, 3, 4, 5, 9 
-1,  otherwise 


(7.9) 


and  the  sequence  of  length  11  is 


+ + - + + + + - . 

Except  for  m • 1 or  2,  Perron  sequences  are  not  m-sequences.  The 
Perron  sequence,  like  any  other  periodic  sequence,  can  be  generated 
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by  some  shift  register  generator,  possibly  with  a nonlinear  feedback 
logic.  Other  binary  pulse  compression  codes  have  been  proposed, 
based  on  dividing  N into  residue  classes  [3,  8].  These  codes  have 
been  generated  for  code  lengths  which  are  prime  numbers,  and  Table 
7.3  shows  the  lowest  achievable  peak  side  lobe  levels  using  these 
codes.  It  has  been  claimed  that,  for  a given  code  length,  such  a 
code  produces  the  lowest  possible  side  lobe  level. 

TABLE  7.3  SIDE  LOBE  LEVELS  FOR  RESIDUE  CLASS  SEQUENCES  [8] 


Code 

Length 

Lowest 
Achievable 
Side  Lobe 

Code 

Length 

Lowest 
Achievable 
Side  Lobe 

Code 

Length 

Lowest 
Achievable 
Side  Lobe 

7 

1 

41 

4 

79 

6 

11 

1 

43 

4 

83 

6 

13 

1 

47 

4 

89 

6 

17 

2 

53 

5 

97 

7 

19 

2 

59 

5 

101 

6 

23 

3 

61 

5 

103 

8 

29 

A 

67 

5 

107 

7 

31 

3 

71 

5 

109 

8 

37 

4 

73 

6 

113 

7 

8.  Complementary  Sequences 


A set  of  complementary  sequences  [11,  23,  60]  is  defined  as  a pair 
of  equally  long,  finite  sequences  of  two  kinds  of  elements  which  have 
the  property  that  the  number  of  pairs  of  like  elements  with  any  given 
separation  in  one  series  is  equal  to  the  number  of  pairs  of  unlike 
elements  with  the  same  separation  in  the  other  series.  When  the  two  ele- 
ments of  the  series  are  taken  to  be  +1  and  -1,  it  follows  immediately 
that  the  sum  of  their  two  respective  autocorrelations  is  zero  except 
at  zero  lag.  For  example,  the  two  sequences 

A»-++-+-+++-  (8.1) 

and  B = -+  + + + + + (8.2) 

are  complementary  since  A has  3 pairs  of  like  adjacent  elements  and  B 
has  3 pairs  of  unlike  adjacent  elements,  A has  4 pairs  of  like  alter- 
nate elements  and  B has  4 pairs  of  unlike  alternate  elements,  and  so 
on.  The  autocorrelations  of  A and  B can  be  easily  seen  to  be 

= (1,  -2,  -1,  2,  1,  0,  -1,  0,  -3,  10,  -3,  0,  -1,  0,  1,  2,  -1,  -2,  1) 

(8.3) 

Pg  = (-1,  2,  1,  -2,  -1,  0,  1,  0,  3,  10,  3,  0,  1,  0,  -1,  -2,  1,  2,  -1) 

(8.4) 

and 

P^  + Pg  = (0,  0,  0,  0,  0,  0,  0,  0,  0,  20,  0,  0,  0,  0,  0,  0,  0,  0,  0). 

(8.5) 

The  necessary  conditions  for  the  existence  of  complementary  sequences 
are  [23]:  (1)  that  the  code  length  N be  even,  and  (ii)  that  N must  be 

expressible  as  a sum  of,  at  most,  two  squares.  Existence  of  complementary 
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sequences  has  been  established  when  N Is  a power  of  two.  Otherwise » 
they  have  been  discovered  only  for  N of  the  type 

N - 10  • 2“  , for  n ^ 0 . 

The  following  are  the  general  properties  of  complementary  sequences: 

(1)  The  numbers  of  elements  In  the  two  sequences  are  equal. 

(11)  The  two  complementary  sequences  are  Interchangeable. 

(Ill)  The  order  of  elements  of  either  or  both  the  sequences  may 
be  reversed  without  affecting  the  complementary  property. 

(Iv)  Each  element  of  one  or  both  the  sequences  may  be  multiplied 
by  -1. 

(v)  Alternate  elements  of  both  the  sequences  may  be  multiplied 
by  -1  without  affecting  the  complementary  property. 

Thus  a single  pair  of  complementary  sequences  can  be  the  basis  for  the 
construction  of  several  pairs  of  complementary  series  using  the  above 
properties. 

Beginning  with  a few  pairs  of  complementary  sequences.  It  Is  pos- 
sible to  construct  a variety  of  longer  series  by  simple  combinations 
[23,  42,  62].  Let  A - (Sj , a^ , ••••,  a^^)  and  B » (bj , b^ , ••••,  b^j) 
be  complementary  with  all  their  elements  either  -fl  or  -1.  Also,  let 

b^  ■ - bj  and  B'  denote  the  series  obtained  by  multiplying  every  ele- 

^k,  ^ 

ment  by  -1.  In  addition.  If  Cj^  Is  -1,  then  A Is  taken  to  be  A and 
A^^  ■ A If  Cj^  ■ +1.  Longer  complementary  sequences  can  then  be  formed 
as  follows: 

(1)  The  two  series  Sj  and  Sj  of  length  2N,  where 
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S^  = (A,B)  = (a,,  a^,  •••*,  , b, , , 


• "n> 


(8.6) 


and 


Sj  « (A,B*)  “ (a,,  a^,  •••',  a^^  , bj,  b^ , 
are  complementary. 


' *>N> 


(8.7) 


(ii)  The  Interleaved  series  and  Tj  of  length  2N,  where 
T,  = (a,,  b,,  a^,  b^,  , a^^,  bj^) 

and 


‘»N> 

(8.8) 

(8.9) 

are  complementary. 


(ill)  If  C = (cj,  Cj,  , Cj^)  and  D * (dj,  dj,  , d^)  are 

complementary  (note  N M is  not  required) , then  the  two 
two  series  Uj  and  of  length  2NM  irfiere 

Cl  c,  Cu  d,  d, 

U,  = (A  *,  A % ••••,  A ” , B \ B % 

and 


, B **) 


(8.10) 


c!  cl 

Uj  = (A  ",  ••••,  A ^ A *,  B " , B B ‘) 


(8.11) 


are  complementary.  Similarly,  the  two  series 


c.  d,  c„  d,  dw 

V,  = (A  ‘ , B ‘ , A ^ B S • • • • , A ”,  B *^) 


(8.12) 


and 


(A  ”,  B ”,  A"^  B ^ A *,  B *) 


1 


(8.13) 


are  also  complementary. 


(Iv)  Let  N be  of  the  type  N ■ 2".  Let  the  two  sequences  be 

.1 


partitioned  Into  L 2 sequences  (S,  i n)  of  equal  length 
such  that 


J 
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A - (A,,  Aj,  , A^)  (8.14) 

and 

B - (Bj,  Bj,  •••.,  Bj^)  . (8.15) 

Then  the  two  series  Wj  and  W2  of  length  2N,  where 

Wj  - (A,,  Bj,  Aj,  Bj,  ••••,  A^,  B^)  (8.16) 

and 

Wj  = (A,,  b;,  Aj,  b;,  •••',  Aj^,  B^)  (8.17) 

are  complementary.  Note  that  H ° 0 yields  the  series  Sj 
and  Sj  whereas  8,  = n produces  the  series  Tj  and  Tj  . 

Thus  using  the  above  construction  procedure  and  the  general  proper- 
ties listed  earlier,  a variety  of  complementary  sequences  can  be  generated. 
It  must  be  noted,  however,  that  even  though  the  autocorrelations  of  two 
complementary  series  add  up  to  zero  at  non-zero  delays,  the  Individual 
autocorrelation  structures  and  the  cross-correlation  between  the  two 
series  vary  widely  between  different  series  of  the  same  length.  Four 
methods  have  been  suggested  for  the  utilization  of  complementary  sequences 
in  pulsed  radars. 

Method  1 [191 ; Here  the  two  complementary  sequences  are  transmitted 
simultaneously.  The  two  received  waveforms  are  then  correlated  using 
two  different  matched  filters,  each  matched  to  one  of  the  sequences. 

The  outputs  of  the  two  filters  are  then  added  to  yield  the  final  out- 
putput.  In  the  absence  of  noise  and  doppler  mismatch,  the  outputs  of 
the  two  matched  filters  are  the  autocorrelations  of  the  two  sequences, 
and  hence  the  final  output  will  have  one  peak  at  a delay  corresponding 


to  the  target  range  and  is  zero  elsewhere. 
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Method  2 [421;  Instead  of  simultaneous  transmission,  the  two  sequences 
of  a complementary  pair  may  be  alternated  on  successive  PRF  Intervals. 
The  receiver  contains  two  matched  filters  matched  to  each  of  the  se- 
quences, and  the  received  signal  is  passed  through  the  filter  corre- 
sponding to  the  code  transmitted  in  that  interval.  The  matched  filter 
outputs  for  several  PRF  intervals  are  then  added  together  (Integrated) . 
In  the  absence  of  noise  and  target  motion,  assuming  that  an  even  number 
of  pulses  are  integrated,  the  output  contains  a single  peak  and  Is  zero 
elsewhere.  However,  the  cancellation  performance  is  extremely  doppler 
sensitive  and  doppler  correction  must  be  applied.  Also,  the  difference 
in  target  range  between  two  Intervals  results  in  a relative  shift  at 
the  matched  filter  outputs,  contributing  to  poor  side  lobe  cancellation. 
Thus  It  is  necessary  to  compensate  for  the  target  motion  both  in  range 
and  in  doppler.  If  the  integration  Is  performed  by  an  analog  device, 
its  frequency  response  is  also  Important.  Rather  than  alternating  the 
two  sequences,  equal  numbers  of  the  two  codes  may  be  used  in  different 
sequencing  patterns.  Different  patterns  alter  the  frequency  spectrum 
of  the  side  lobe  energy  making  it  more  or  less  suitable  for  the  par- 
ticular device  being  used  for  integration.  It  is  also  conceivable  that 
several  different  pairs  of  complementary  sequences  could  be  used  within 
an  integration  period. 

The  two  methods  described  above  were  proposed  for  a general  radar 
application.  The  following  two  methods  are  designed  specifically  for 
tracking  applications  where  large  side  lobes  are  tolerable  if  they  are 
far  away  from  the  main  lobe. 


Method  3 [111;  Let 

A - (Sj , aj,  , a^^) 


(8.18) 


and  B » (bj,  bj,  **",  bj^)  (8.19) 

be  a pair  of  complementary  sequences.  Consider  the  sequence  C of 
length  2N+K  such  that 
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C » (Cj,  Cj, 


‘^2N+K^ 


(8.20) 


where 

Cj  = for  1 < N 

Cj  * 0 for  N <ii  N+K  (8.21) 

and  c,  * b,  „ „ for  N+K  < 1 ^ 2N+K 
I 1— N-K 

Using  the  fact  that  A and  B are  complementary,  the  autocorrelation  of 
C can  be  easily  shown  to  be 


P (0)  = 2N 
c 

P^(j)  “ 0 for  1 < j < K 

P (j)  - P.  (N-j+K)  for  K ^ j ^ 2N+K 
c oa 

where  P,  denotes  the  cross-correlation  between  B and  A. 
ba 


(8.22) 


For  example. 


let  A * (1,  1,  -1,  1)  and  B = (1,  1,  1,-1),  N = 4 and  K = 5.  Then, 

C ■=  (1,  1,  -1,  1,  0,  0,  0,  0,  0,  1,  1,  1,  -1)  (8.23) 

and  - (-1,  0,  3,  0,  1,  0,  1,  0,  0,  0,  0,  0,  8,  0,  0,  0,  0,  0,  1,  0,  1,  0,  3,  0,  -1)  . 

(8.24) 

Thus  the  effect  of  Inserting  5 zeros  between  A and  B is  to  achieve  an 
autocorrelation  with  zero  side  lobes  in  the  interval  [-5,  5].  Figure 
8.1(a)  and  8.1(b)  depict  the  signal  and  the  corresponding  autocorrela- 
tion function.  As  can  be  seen,  the  autocorrelation  is  precisely  the 
type  desired  for  cracking  applications.  The  main  disadvantage  of  such 


(b)  Aut 

Effect  of 
Series 
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I 

1 


1 

I 
j 

a waveform  Is  that  the  signal  power  utilization  efficiency  decreases  I 

as  K increases  since  no  signal  is  transmitted  during  the  period  between 

the  two  pulses.  However,  if  the  radar  is  an  electronically-steered 

multipurpose  radar,  it  is  possible  for  the  radar  to  perform  some  other 

function  during  the  period  between  the  pulses.  In  such  an  application, 

by  selecting  K large  enough,  the  signal  may  also  be  useful  during 

search  operations. 

Method  4 [111;  Let  A and  B be  a pair  of  complementary  sequences  of 

length  N and  consider  the  code  C of  length  4MN  formed  as  follows: 

C = (A,  B,  A',  B,  ) . (8.25) 

repeated 
M times 

The  peak  value  of  the  autocorrelation  occurs  at  zero  lag  and  is  given  by 

P^(0)  - 4MN  . (8.26) 

Since  the  code  is  formed  by  a periodic  repetition,  P^(j)  must  have 
large  side  lubes  when  J is  a multiple  of  the  period  4N.  Using  the 
property  that  A and  B are  complementary,  the  autocorrelation  of  C in 
the  region  0 < j < 4M  can  be  written  as  follows: 

Pc<J)  “ - Pgb^”  “ for  0 < j S N 

PcO)  - - Pab<N- j)  + (2M-1)  [P^(2N- j)-Pg(2N- j)]  for  N < j S 2N 

P^.(J)  - + Pab<J  - 3N)  + (2M- 1)  tP|,(2N- j) -P^(2N- j)]  for  2N  < j S 3N 

Pc(j)  - Pab<J  - 3N)  for  3N  < j i 4N  . (8.27)  j 


Therefore,  in  the  region  0 < j i N,  the  autocorrelation  of  C is  governed 
by  the  cross-correlation  between  the  sequences  A and  B.  For  the  tracking 
application,  let  the  region  of  interest  be  -K  i j i K where  K i N,  and  let 
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q 


ba 


» nsx 

Jik 


(N  - J)1  ; 


(8.28) 


then  the  peak-to-side  lobe  ratio  in  the  region  [-K,  K]  is 


P (0) 
c 


^ba 


4t1N 

%a 


(8.29) 


This  quantity  has  a lower  bound  of  4M  and  can  be  maximized  by  minimizing 

q . Therefore,  the  best  pair  of  complementary  sequences  is  one  where 
D3 

the  cross-correlation  is  small  in  the  region  of  interest.  One  must  note 
here  that  the  region  of  low  side  lobes  cannot  be  extended  beyond  j = N 
because  of  the  autocorrelation  structure.  Other  series  combinations  of 
complementary  sequences  could  be  considered  but  no  useful  results  have 


been  published. 
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9.  Simultaneous  Amplitude  and  Phase  Modulation 
(13,  28,  33,  56] 


In  all  of  the  codes  considered  In  the  previous  sections,  the  pulse 
amplitude  was  held  constant  while  its  phase  was  varied.  The  Huffman 
codes  [11,  33]  are  made  up  of  finite  sequences  of  amplitude  and  phase 
modulated  pulses.  The  complex  representation  of  such  a signal  can  be 
written  for  an  arbitrary  N as 


I|»(t)  =•  u(t)  e-^ 

N 

where  u(t)  = a^R^(t) 

n = 1 


(9.1) 


(9.2) 


R (t)  = 1 for  (n  - 6)  < t < n6 
n 


= 0 elsewhere 


and 


N6  . 


(9.3) 


(9.4) 


The  excitation  sequence  (Sj,  a^,  a^^)  consists  of  complex  numbers 

whose  magnitudes  and  phases  are  not  necessarily  equal.  The  magnitude 
variation  provides  the  amplitude  modulation  and  the  phase  variation 
results  In  the  phase  modulation.  The  elements  of  the  complex  excita- 
tion sequence  are  chosen  In  such  a way  that  the  autocorrelation  has  a 
desirable  structure.  In  the  case  of  Huffman  codes,  the  autocorrelation 
of  the  sequence  A • (a^,  a^,  *•*•,  3j^)  is  restricted  to  be 

■E.  2 = 

k=  1 

= 0 for  0 < |j|  < N - 1 (9.5) 


P3(0) 


P (j) 
a 


P (N  - 1) 
a 


± 1 
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and  P ( - (N-  D)  = a, a*  = ± 1 
a 1 N 

where  E Is  Che  energy.  Because  of  this  autocorrelation  structure,  these 

sequences  are  also  called  Impulse  equivalent  sequences.  If  A(z)  Is  the 

Z-transform  of  {a  1 and  P(z)  Is  the  Z-transform  of  the  autocorrelation 
n 

P (j),  then  It  can  be  easily  shown  that 

iX 

P(z)  = A(z)  A*(z  *)  (9.6) 


where 


A(z)  » 

+ a^z  ^ 4-  ajZ 

2 

+ • • 

, -(N- 

+ 

■1) 

(9.7) 

and 

N-1 

P(z)  = 

-.1 

(9.8) 

-(N-1) 

Substituting 

, for  P (j)  from 

iX 

(9.5), 

P(z)  becomes 

P(z)  - 

a,  a*  z*^"^  + E + 

(N- 

-1) 

(9.9) 

The 

2N-2  roots 

of  (9.9)  can 

be  shown 

Co  be  as 

follows : 

(1) 

If 

1 « 

for  0 < 

k < 

N - 2 

(9.10) 

-N+1)(J) 

for  N - 

1 < 

k < 2N-2 

(9.11) 

(11)  If  a*aj^  = + 1 , 


-1 


j^^^j(k+0.5)<^ 

_ J(k-N+1-H).5)4. 

A w 


for  0 ^ k < N - 2 
for  N - 1 < k S 2N  - 2 


where 

♦ 


2tt 


(9.12) 

(9.13) 


N-  1 


(9.14) 
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H»-‘  - f * 


" R, 


^ - 1 
2 


(9.15) 

(9.16) 


Therefore,  the  roots  lie  on  two  concentric  circles  centered  at  the 

origin  and  are  equally  spaced  around  the  circles.  There  are  N-1  pairs 

of  roots  where  each  pair  consists  of  one  root  on  the  inner  circle  and 

one  on  the  outer  circle  along  the  same  radial  line.  One  of  each  pair 

of  these  roots  can  be  assigned  to  A(z)  while  the  other  is  automatically 

assigned  to  A*(z  *).  Once  the  roots  of  A(z)  are  known,  the  corresponding 

polynomial  and  hence  the  excitation  sequence  can  be  determined.  Since 

N-1 

one  from  each  of  N-1  pairs  are  assigned  to  A(z),  there  are  2 dif- 
ferent possible  choices  of  these  roots.  While  every  one  of  these 
choices  results  in  the  same  autocorrelation,  some  yield  real  sequences 
and  the  others  result  in  complex  sequences.  In  addition,  each  choice  re- 
sults in  a different  distribution  of  the  amplitudes.  For  maximum  trans- 
mission efficiency,  the  radar  should  transmit  at  peak  power  over  the 
entire  duration;  i.e.,  the  ideal  situation  is 


laj  = UJ  la^jl  = • (9.17) 

The  choice  of  the  root  pattern  must  be  such  that  the  amplitude  distri- 
bution is  as  close  as  possible  to  the  ideal  distribution  (9.17)  . That 
is,  the  root  pattern  must  be  chosen  to  maximize  the  energy  utilization 
ratio  [1,  2,  50], 


(9.18) 
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A computer  evaluation  has  been  carried  out  to  determine  the  optimum  root 
pattern  [36].  It  has  been  shown,  for  example,  that  for  N = 11,  the 
highest  possible  value  of  a Is  0.54. 

In  order  to  Improve  the  energy  utilization  ratio,  the  restriction 
on  the  autocorrelation  structure  can  be  relaxed.  For  example,  if  an 
additional  side  lobe  pair  of  amplitude  one  is  allowed  to  exist  one  time 
side  lobe  width  away  from  the  unavoidable  end  time  side  lobe  pair  [35], 
i.e. , 

P^(j)  = E for  j = 0 

a 

iPa^^^l  = 1 for  |j|  = N - 1 

(9.19^ 

|P  (j)|  = 1 for  111  = N - 3 

a 

P (j)  =0,  otherwise 

then,  for  N * 11,  the  energy  utilization  ratio  can  be  increased  to  0.624. 
Thus,  in  tracking  applications,  by  allowing  even  more  side  lobes  to  exist 
outside  the  region  of  interest,  higher  energy  utilization  ratios  could  be 
achieved. 

Another  approximation  to  the  Huffman  code  is  the  modified  Barker 
code  [34]  where  the  autocorrelation  is  specified  as 

P^^(j)  = E for  j * 0 

lP^(j)|  = 1 for  |j|  = N - 1 

(9.20) 

l^g(j) I = P for  |j j » 2,4,6,  

P (j)  =0,  otherwise 

Q. 

Clearly,  p * 0 yields  the  Huffman  code  and  p = 1 gives  the  Barker  code 
If  it  exists  for  the  given  length.  The  code  can  be  generated  for 


different  values  of  p using  a similar  procedure  as  the  Huffman  code. 

The  energy  utilization  ratio  Increases  as  p increases  from  0 to  1.  For 
example,  if  N ••  13,  p • 0 gives  a value  of  o * 0.564  whereas  p = 0.5 
yields  o = 0.713. 

Ternary  sequences  are  amplitude  and  phase  modulated  codes  where 
both  the  amplitude  and  phase  take  on  one  of  two  values.  Specifically, 
a ternary  sequence  {a^}  consists  of  elements  1,  -1,  and  zero.  Ternary 
sequences  whose  cyclic  autocorrelation  is  zero  everywhere  except  at 
zero  lag  have  been  considered.  Tompkins  [56]  managed  by  brute  force 
technique  to  generate  codes  up  to  length  18.  A systematic  method  for 
generating  ternary  sequences  is  described  by  Chang  [13]  when  the  length 
of  the  sequence  Is  given  by  N = (3"  - 1)  where  n is  odd.  The  non- 
periodic autocorrelation  of  these  functions  and  the  effect  of  doppler 
mismatch  have  not  been  investigated. 
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10.  Binary  Coding  of  In-Phase  and  Quadrature  Components 


The  complex  envelope  defined  in  (3.3)  may  be  written  as 
u(t)  - a(t) 


X (t)  - Jx  (t) 
d q 


(10.1) 


where  x^(t)  and  *^(1)  are  the  In-phase  and  quadrature  components, 
respectively.  The  corresponding  transmitted  signal  is 


s(t;o)g)  = x^(t)  cosWgt  + X (t)sln(jgt 


(10.2) 


If  both  Xjj(t)  and  x^(t)  are  binary  coded,  l.e. , if  they  take  on  values 
+1  and  -1,  then  it  can  be  easily  seen  that 

a(t)  * ^2 


and  (Ht)  = 1 ^ 3 4 ’ 


(10.3) 

(10.4) 


which  represents  a four-phase  code  of  constant  amplitude.  Conversely, 
any  four-phase  code  where  the  phases  are  equally  separated  can  be  im- 
plemented by  transmitting  a pair  of  binary  coded  signals  in  quadrature. 
Since  good  four-phase  codes  of  length  greater  than  16  are  not  known, 

various  combinations  of  binary  codes  must  be  tried  for  x,(t)  and  x (t) . 

d q 

The  autocorrelation  of  u(t)  can  be  expressed  in  terms  of  the  auto- 
and  cross-correlation  Integrals  of  the  in-phase  and  quadrature  components; 
l.e. , 

Puu<^>  “ ^Pdd^'"^  Pqq^">>  - j^Pdq^^>  " Pqd^''^^  * 

The  radar  performance  depends  upon  the  peak-to-slde  lobe  ratio  of  |p  (t) | , 


uu 


uu 


qq 


dq' 


J a ' 

qd 


The  problem  Is  to  select  the  binary  sequences  and  x^  which  minimize 

for  6 i T i k6  where  k represents  the  window  of  Interest.  At 

the  present  time  there  is  no  rational  method  for  determining  x^  and  x^. 

Many  pairs  of  (x.,  x ) were  considered.  Good  results  were  obtained  for 
d q 

symmetric  pairs,  in  which  case  P^q('r)  " ~ example, 

X.  - + + + + — + + -+-  + -++-  + + -+-+-++--+  + + + 

a 

X * ++ + + --  + - + -+--++ + + 

q 

Ip^^I  = (62,  0,  2,  0,  2,  0,  2,  0,  6,  •••)  . 

More  work  is  required  before  this  technique  can  be  evaluated. 


11.  Example  of  Binary  Code  Selection 


As  explained  earlier,  the  erlterlon  for  the  selection  of  a code 
depends  upon  the  application.  The  PSLR  or  the  peak  - to  - rms  side  lobe 
ratio  is  appropriate  for  a search  radar.  On  the  other  hand,  for  a radar 
which  is  dedicated  to  acquisition  and  tracking,  only  the  side  lobes  In 
tiie  vicinity  of  the  main  lobe  are  of  concern.  For  example,  the  selection 
of  the  signal  may  be  based  upon  its  autocorrelation  at  lags  0,  6,  26,  36, 
46,  and  56  which  corresponds  to  a total  of  10  range  cells  ( one  range  cell 
is  equivalent  to  the  width  of  the  compressed  pulse,  6).  The  following  Is 
a comparison  of  several  binary  codes  of  length  31,  where  PPO  stands  for 
peak-to-side  lobe  ratio  In  dB  for  the  entire  correlation,  PRO  is  the  over- 


all peak-to-rms  side  lobe  ratio  in  dB,  PPL  and  PRL  are  the  corresponding 
quantities  witiiin  5 range  cells  of  tlie  main  peak. 

(1)  Contiguous  Barker  Codes  : A code  of  length  31  may  be  formed 

by  arranging  Barker  codes  of  length  13,  13,  and  5 contiguously. 

Tlie  performance  figures  for  such  a code  are  PPO  = 5.7,  PRO  = 18.2 
PPL  = 20.3,  and  PRL  = 21.5.  Contiguous  arrangement  of  Barker 
codes  of  lengths  13,  11,  and  7 also  yields  a code  of  length  31. 
For  this  sequence,  PPO  = 9,  PRO  = 18.1,  PPL  = 20.3,  PRL  = 26.8. 

(11)  Linear  m - sequences  : There  are  3 feedback  connections  and  31 

possible  Initial  conditions  for  each  connection,  yielding  a 
total  of  93  linear  m - sequences.  Only  those  whose  performance 
is  optimum  in  some  sense  are  mentioned  here.  There  are  17 
sequences  with  the  highest  PPO  of  17.8,  the  feedback  connections 


and  initial  conditions  of  which  are  given  in  table  7.2.  The 
sequence  with  feedback  15,4,3,2)  and  Initial  condition  (0,1, 0,0,1) 
has  the  highest  PRO  wlilch  is  PRO  = 25.  The  sequence  with  the 
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highest  PPL  and  E'RL  is  obtained  by  the  feedback  connection 
[5,2]  and  Initial  condition  (0,1, 1,0,0),  and  the  values  are 
PPL  = 29.8  and  PRL  = 33.8. 

(Ill)  Residue  Class  Sequences  : Of  the  many  residue  class  sequences 
available  for  N = 31,  only  a few  were  analyzed.  The  important 
result  Is  that  a sequence  generated  by  dividing  31  into  10 
residue  classes  gives  a PPO  of  20.3  wliich  is  larger  than  that 
of  tile  best  m - sequence  and.  In  fact,  is  the  liighest  that  can 
be  obtained  for  a binary  code  of  length  31.  None  of  the  other 
sequences  studied  li.id  larger  values  for  PPL  than  the  best 
m - sequence  given  above. 

(tv)  Perron  Sequences  : The  residue  class  sequence  obtained  using 
two  residue  classes  is  identical  to  the  Perron  sequence.  No 
such  sequence  was  found  which  had  better  performance  than  the 
m - sequences. 

(v)  Larger  Length  Sequences  ; A sequence  of  lengtli  35  can  be 

obtained  by  combining  barker  codes  of  lengths  7 and  5.  Such 
a sequence  has  PPI,  = 16.9  and  PRL  = 19.8.  A residue  class 
sequence  of  length  N = 47  was  also  studied.  This  sequence 
has  the  highest  PPL  .and  PRL  of  all  the  sequences  considered 
and  the  values  are  PPL  = 33.4  and  PRL  = 37.4. 

Table  11.1  shows  those  sequences  which  are  optimum  in  some  sense. 

The  performance  of  these  sequences  may  be  improved  further  by  the  use 
of  a weighting  filter. Such  a filter  Is  especially  simple  to  synthesize 
for  the  combined  Barker  code.  The  behavior  of  the  binary  sequences  in 
tlie  presence  of  doppler  mismatch  w.is  not  investigated  in  detail. 


12.  ■Suniro.try  and  Conclusions 


The  following  types  of  pulse  compression  signals  were  discussed  in 
this  report:  (i)  linear  and  nonlinear  FM,  (ii)  polyphase  codes,  (iii) 

Barker  codes,  (iv)  shift  register  generated  binary  sequences,  (v)  com- 
plementary series,  and  (vi)  codes  with  simultaneous  amplitude  and  phase 
modulation.  The  performance  of  a pulsed  radar  is  determined  by  exami- 
nation of  the  ambiguity  function  when  a single  pulse  is  involved.  The 
performance  may  be  improved  either  by  integration  or  by  weighting  or 
by  both.  In  the  evaluation  of  a pulse  compression  signal,  the  follow- 
ing factors  will  be  considered:  (1)  signal  generation,  (ii)  peak-to- 

slde  lobe  ratio  (PSLR)  of  the  autocorrelation  and  the  relative  location 
of  the  first  large  side  lobe,  (ill)  effect  of  doppler  mismatch,  and 
(iv)  applicability  of  integration  and/or  weighting. 

Linear  FM  signals  can  be  generated  for  BT  products  up  to  4000. 

The  PSLR  is  about  13.5  dB  which  implies  that  weighting  is  required. 

The  first  side  lobe  is  the  largest.  A doppler  mismatch  causes  both  a 
reduction  in  the  PSLR  and  an  error  in  range  measurement.  Using  weight- 
ing techniques,  the  PSLR  can  be  increased  at  the  expense  of  a small 
widening  of  the  main  lobe.  For  example,  Hamming  weighting  can  increase 
the  PSLR  to  42.8  dB  with  a main  lobe  widening  factor  of  1.47.  The  PSLR 
can  be  increased  by  using  nonlinear  FM  but  the  performance  becomes  very 
sensitive  to  doppler  mismatch. 

Frank  polyphase  codes  are  the  most  important  polyphase  codes  and 
exist  for  lengths  N = p^ . For  large  values  of  p,  the  PSLR  approaches 
Tip.  For  example,  p * 6(N  = 36)  corresponds  to  PSLR  = 25  dB.  The 
polyphase  codes  have  doppler  degradations  somewhat  similar  to  those  of 


linear  FM.  Another  interesting  type  of  polyphase  code  is  the  generalized 
Barker  sequence  where  the  PSLR  = N,  the  length  of  the  code.  Unfortunately 
the  existence  has  not  been  established  for  lengths  larger  than  15. 
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The  most  useful  binary  codes  are  Barker  sequences,  shift  register 
generated  sequences,  and  complementary  series.  Barker  codes  exist  only 
for  N = 2,3,4,5,7,11,  and  13  and  have  the  unique  property  that  the  PSLR 
la  equal  to  N.  The  performance  deteriorates  rapidly  for  large  values  of 
doppler  mismatch.  Weighting  filters  can  be  designed  to  Increase  sub- 
stantially the  PSLR.  Particularly  important  to  tracking  applications 
is  the  fact  that  the  first  K side  lobes  can  be  entirely  suppressed  by 
a transversal  filter  of  order  2K.  Longer  codes  can  be  achieved  using 
multistage  Barker  codes. 

Linear  m-sequences  are  commonly  used  because  they  are  easily  gen- 
erated by  linear  shift  register  generators.  An  n-stage  shift  register 
with  proper  feedback  can  generate  an  m-sequence  of  length  N = 2”  - 1. 

Since  the  signal  is  aperiodic,  the  ambiguity  function  depends  upon  the 
feedback  configuration  and  on  the  initial  condition  for  a given  feed- 
back structure.  Initial  conditions  which  produce  the  largest  PSLR  have 
been  tabulated.  For  example,  the  best  m-sequence  of  length  31  has  a 
PSLR  of  12.9  dB.  For  large  values  of  N,  the  PSLR  approaches  /¥  , which 
can  be  further  Increased  by  the  use  of  a weighting  filter  similar  to  the 
one  employed  in  the  case  of  a Barker  code.  The  performance  of  m-sequences 
in  the  presence  of  doppler  mismatch  is  superior  to  that  of  the  Barker 
code  since  it  has  no  appreciable  spurious  peaks. 

Another  type  of  binary  sequence  is  the  Perron  sequence  whose  periodic 


correlation  is  maximum  at  zero  delay  and  is  minus-one  elsewhere.  Such 


sequences  exist  for  lengths  of  the  type  N * 4m  - 1 when  this  number  is 
prime.  Unfortunately,  there  Is  no  general  theory  about  the  aperiodic 
correlation  properties  of  the  Perron  sequences. 

Residue  class  sequences  have  been  generated  for  code  lengths  which 
are  prime  numbers.  It  has  been  claimed  that  such  codes  produce  the 
largest  possible  PSLR  for  a given  code  length.  A definite  advantage  of 
these  sequences  over  the  m-sequences  Is  that  they  exist  for  a larger 
set  of  lengths. 

It  Is  possible  to  obtain  longer  binary  sequences  by  combining  two 
or  more  binary  sequences.  It  is  advantageous  to  start  from  a pair  of 
complementary  sequences.  For  example,  a code  can  be  formed  by  repeat- 
ing M times  (A,  B,  A',  B)  where  A and  B are  complementary  of  length  N. 
The  PSLR  within  K (K  < N)  steps  of  the  main  lobe  has  a lower  bound  of 
4M  and  can  be  maximized  by  minimizing  the  cross-correlation  between 
A and  B.  Another  proposal  was  to  form  a long  code  by  Inserting  K zeros 
between  two  complementary  sequences  of  length  N.  In  this  case,  the 
side  lobes  will  be  zero  within  K steps  of  the  main  lobe,  which  Is  par- 
ticularly attractive  In  tracking  applications. 

Improved  performance  may  be  obtained  by  Integration  of  two  or 
more  output  pulses.  If  the  matched  filter  outputs  due  to  two  comple- 
mentary sequences  are  Integrated,  the  result  will  have  a single  peak 
and  will  be  zero  elsewhere  (assuming  zero  doppler  mismatch) . Two 
methods  were  suggested  to  take  advantage  of  this  property.  In  the 
first  method,  the  two  complementary  sequences  are  transmitted  simul- 
taneously. The  two  received  waveforms  are  then  passed  through  the 
matched  filters,  each  matched  to  one  of  the  sequences,  and  their  outputs 
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are  added  Co  yield  Che  final  ouCpuC.  In  Che  second  mechod,  Che  Cwo 
sequences  of  a complemenCary  pair  are  alCernaCed  on  successive  PRF  InCer- 
vals  and  Che  maCched  fllCer  ouCpuCs  for  several  Incervals  are  inCegraCed. 
In  chls  case,  Che  performance  is  exCremely  doppler  senslclve  and  doppler 
correcCion  musC  be  applied. 

In  all  Che  signals  considered  Chus  far,  higher  BT  produces  were 
accomplished  by  phase  modulaClon.  An  alCernaCive  scheme  Co  achieve  high 
BT  produce  Is  slmulCaneous  ampllCude  and  phase  modulaClon.  An  example 
of  such  a signal  Is  Che  Huffman  code  which  exlsCs  for  all  lengths.  If 
Che  lengch  of  Che  code  is  N,  Chen  Che  auCocorrelaClon  is  maximum  ac  zero 
lag  and  Is  zer  ) elsewhere  excepC  ac  lag  N-1.  Several  Huffman  sequences 
are  available  for  a given  lengCh,  buC  Che  besC  sequence  is  one  which 
resulCs  in  Che  maximum  energy  uCllizaClon  racio,  i.e.,  which  has  Che 
mosC  uniform  ampllCude  disCribuCion.  In  cracking  appllcaClons,  oCher 
side  lobes  may  be  allowed  Co  exlsC  In  addlClon  Co  Che  ones  In  the 
Huffman  code.  Such  sequences  can  be  generated  by  a procedure  similar 
to  the  one  used  for  Huffman  codes  and  result  in  higher  energy  utiliza- 
tion ratios  as  compared  to  the  Huffman  codes.  For  example,  for  N = 11, 
the  highest  energy  uClllzatlon  ratio  for  a Huffman  code  is  0.54,  whereas 
by  allowing  one  more  side  lobe  at  lag  N - 3 - 8,  the  ratio  can  be  In- 
creased to  0.624. 

The  emphasis  of  Chls  study  has  been  on  binary  phase  coded  signals 
suitable  for  Cracking  appllcaClons.  As  mentioned  earlier,  a signal  can 
be  used  In  Cracking  radars  even  if  its  autocorrelation  has  a low  PSLR 
as  long  as  the  side  lobe  levels  are  small  in  Che  vicinity  of  Che  main 
lobe.  In  view  of  this  criterion,  Che  following  conclusions  are  drawn 
regarding  binary  phase  coded  signals: 


74 


(1)  Long  codes  and  hence  high  compression  ratios  can  be  obtained 
by  combining  Barker  codes.  The  advantage  of  these  codes  Is 
that  weighting  filters  can  be  easily  designed  to  suppress 
the  side  lobes  near  the  main  lobe.  The  critical  considera- 
tion In  evaluating  such  a code  Is  Its  sensitivity  to  doppler 
mismatch. 

(11)  Linear  m-sequences  are  easy  to  generate  and  are  less  sensi- 
tive to  doppler  shift  than  Barker  codes.  Based  either  on 
theoretical  considerations  or  on  exhaustive  computer  evalua- 
tion optimum  m-sequences  have  been  tabulated  which  yield 
the  highest  PSLR.  These  codes  must  therefore  be  reevaluated 
for  tracking  applications  where  a high  PSLR  Is  not  essential. 
Furthermore,  weighting  filters  can  be  designed  to  Improve 
their  performance. 

(Ill)  Residue  class  sequences  exist  for  many  more  lengths  than 

the  m-sequences.  Even  though  It  Is  claimed  that  these  se- 
sequences  have  the  best  PSLR  for  a given  length,  their 
suitability  In  tracking  has  not  been  established. 

(iv)  Zero  side  lobes  In  the  vicinity  of  the  main  lobe  can  be 
attained  by  forming  a long  code  where  zeros  are  Inserted 
between  a pair  of  complementary  sequences.  The  main  draw- 
back to  such  a scheme  Is  the  reduced  energy  utilization 
ratio. 

(v)  Long  codes  may  be  obtained  by  repeating  a pair  of  comple- 
mentary series.  In  such  a situation,  the  side  lobes  level 
near  the  main  lobe  Is  determined  by  the  cross-correlation 


between  the  two  complementary  series.  The  optimum  choice 
of  the  complementary  pair  is  therefore  one  that  has  the 
smallest  cross-correlation.  No  theoretical  tool  is  avail- 
able for  this  optimization,  and  hence  it  must  be  accom- 
plished via  a computer  search. 

(vi)  Complementary  series  may  be  used  in  conjunction  with  pulse 
integration  to  yield  an  Improved  performance.  A pair  of 
complementary  series  may  either  be  transmitted  simultaneously 
or  be  alternated  on  successive  PRF  intervals,  and  the  outputs 
Integrated.  In  the  absence  of  doppler  shift,  all  the  side 
lobes  are  completely  suppressed. 

(vll)  It  was  shown  that  a four-phase  modulation  can  be  easily  ob- 
tained by  binary  modulation  of  both  in-phase  and  quadrature 
components.  However,  no  rational  method  for  the  selection 
of  the  pair  of  sequences  has  yet  been  developed. 

(vill)  The  selection  of  binary  codes  of  length  31,  35,  and  47  was 
discussed,  resulting  in  the  recommendation  of  optimum  codes 
where  both  local  and  overall  performances  were  considered. 

Of  course,  better  performance  may  be  obtained  by  using  signals  other 
than  binary  phase  coded  signals.  Examples  of  such  signals  are  the  poly- 

I 

phase  codes  and  signals  with  simultaneous  amplitude  and  phase  modulation. 
Invariably  the  Improvement  in  performance  is  attained  at  the  expense  of 
increased  complexity  and  cost  of  implementation. 
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